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Abstract 

We study the harmonic moments of Galton- Watson processes, pos- 
sibly non homogeneous, with positive values. Good estimates of these 
are needed to compute unbiased estimators for non canonical branch- 
ing Markov processes, which occur, for instance, in the modeling of 
the polymerase chain reaction. By convexity, the ratio of the har- 
monic mean to the mean is at most 1. We prove that, for every square 
integrable branching mechanisms, this ratio lies between 1 — A/k and 
1 — A! /k for every initial population of size k > A. The positive con- 
stants A and A' such that A^ A' are explicit and depend only on the 
generation-by-generation branching mechanisms. In particular, we do 
not use the distribution of the limit of the classical martingale asso- 
ciated to the Galton- Watson process. Thus, emphasis is put on non 
asymptotic bounds and on the dependence of the harmonic mean upon 
the size of the initial population. In the Bernoulli case, which is rele- 
vant for the modeling of the polymerase chain reaction, we prove essen- 
tially optimal bounds that are valid for every initial population k 1. 
Finally, in the general case and for large enough initial populations, 
similar techniques yield sharp estimates of the harmonic moments of 
higher degrees. 
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Introduction 



We study the behaviour of the harmonic means 1/E(l/Z n ) of Galton- Watson 
processes (Z n ) n ^o, possibly non homogeneous, with positive values. A moti- 
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vation for this theoretical problem is the construction of unbiased estimators 
for samples of branching Markov processes, when the state of an individual 
depends on the number of its siblings. An instance, outside the realm of 
pure probability, where this construction is needed, arises in the modeling 
of the polymerase chain reaction by branching processes, see Sun (1995). In 
this specific case, the offspring of each individual is 1 or 2, the state of the 
first descendant is identical to the state of its parent and the state of the 
other descendant, if any, is a stochastic function of the state of its parent. 
One wishes to estimate, for instance, the mutation rate of the reaction from 
a uniform sample of a given generation. Any unbiased estimator of the state 
of such a sample requires to compute the harmonic mean size of the cor- 
responding generation. But there exists no closed form of these harmonic 
means, except for small initial populations and for small numbers of genera- 
tions. Since the mean sizes of the generations of a branching process are well 
known, the above problem is usually circumvented by assuming that the ini- 
tial population is very large. Then, an averaging effect occurs which implies, 
roughly speaking, that the harmonic mean size of a generation is close to 
its mean size. In the context of the polymerase chain reaction, we showed 
in previous papers, see Piau (2004), that this approximation is accurate for 
surprisingly small initial populations, and we provided sharp quantitative 
estimates of the discrepancy between the harmonic mean and the mean, for 
any initial population. These results also proved useful to establish rigorous 
confidence intervals for the estimator of the mutation rate of the polymerase 
chain reaction, see Piau (2005). Our purpose in the present paper is to give 
the exact extent of this approximation phenomenon for general, possibly 
non homogeneous, Galton- Watson processes with positive values. When 
the approximation phenomenon indeed occurs, we quantify it through non 
asymptotic and essentially optimal bounds. 



1 Results 

In the following, (Z n ) n ^>o denotes a positive Galton- Watson process, possibly 
non homogeneous. The distribution of this Markov process with values 
in {1,2,...} is characterized by a sequence S := (£ n )n^i of distributions 
on {1,2,...}, as follows. For every n ^ 1, conditionally on the past of 
the process, Z n is the sum of Z n -\ random variables of law £ n which are 
independent of the past. Assume that each £ n is integrable of mean fi n ^ 1. 
Then Z n is integrable and, if denotes the expectation when Zq = k, for 
any positive integer k, 

n 

E k (Z n ) = kM n , withM„ := JJ/ii. 

i=i 
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On the other hand, by convexity, the sequence of general term M n Efc(l/.Z n ) 
is nondecreasing for n ^ 0. Thus every term is at least 1/k. Our aim 
is to provide explicit bounds of the harmonic moments, which imply, in 
particular, that M n K k (l/Z n ) is close to 1/k when this is so. In other words, 
we wish to show that the sequence of general term M„Efc(l/Z n ) is nearly 
constant. Indeed, for every fixed n ^ and when k — > oo, the law of large 
numbers implies that E/%(1/Z n ) is equivalent to l/(k Ei(Z„)), whose value 
is l/Efc(Z n ) = l/(kM n ). Much more is true, as we show below. To ease the 
task of the reader, we first state the consequence of our general results, in 
the homogeneous case. 

Theorem A Assume that E is constant and square integrable. Thus £„, = £ 
and M n = fi n where £ is square integrable and (j, n =: fj, ^ 1 for every n ^ 1. 
Then, there exists a positive constant A, which depends only on £, such that, 
for every integer k > A and every n ^ 0, 

l/k^fi n E k (l/Z n )^l/(k-A). 

Assume furthermore that /j / 1. There exists a positive constant A' , which 
depends only on £ ; such that A' ^ A and, for every integer k > A 1 , 

lim ^ n E k (l/Z n ) ^ l/{k-A'). 

n— too 

1.1 Harmonic moments 

Theorem EI is a consequence of a general quantitative result, stated as the- 
orem El below, which deals with non homogeneous processes. To state and 
prove this result, we rely on some specific families of distributions, that we 
define now. 

Definition 1 For every m ^ 1, the generating function g m of the positive, 
integer valued, random variable L m is such that, for any t E [0, 1], 

E(t Lm ) := g m (t) := t/(m - (m - l)t). 

For any positive c, the random variable L C)Ta is such that, for any t G [0, 1], 

E(t Lc "") := 9c , m (t) ■= (g m (t c )) 1/c = t/(m - (m - l)* c ) 1/c - 

Thus, 

9m — 51, m- For every c, F c ,m ^ 1 almost surely and E(i C)m ) — m. 
When m = 1, L m = L c ^ m = 1 almost surely. When m > 1, the distribution 
of L m — 1 is geometric and the distribution of {L c ,m — l)/c is negative 
binomial. 
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Definition 2 For any positive c, let A c denote the set of distributions of 
integrable random variables L ^ 1 such that, for any t £ [0, 1], 

IE(t L ) < fl c ,m(t), m:=E(L). 

For any positive c, let A' c denote the set of distributions of integrable random 
variables L ^ 1 snc/i i/iai, /or any t E [0, 1], 

E(t L ) ^ g c ,m(t), m := E(L). 

Note that one compares the distribution of L to distributions of random 
variables, not a priori integer valued but with the same mean. We are now 
able to state our main result. 

Theorem B (1) Let n 1. Assume that there exists c such that £j E «4 C 
/or every i ^ n. Then, for every k > c, 

M n E fc (l/Z n ) ^ l/(fc-c). 

(2) Assume that M n — » oo u>/ien n — > oo and f/iaf i/iere exists c such that 
£i G .44 /or et>erj/ i ^ 1 . Then, /or ewery k > c, 

lim M n E fc (l/Z n ) > l/(fc-c). 

n— >oo 

Recall that, by convexity, the sequence M n Efc(l/Z n ) is nondecr easing, hence 
the existence of the limit when n — > oo is a general fact. Assertion (2) 
becomes false when M n is allowed to stay bounded, or when one replaces 
the limit n — > oo by a finite n since, for instance, the n = value is 
On the other hand, in practical situations, the hypothesis that M n — > oo is 
easy to check since it only involves the first moments of the generation-by- 
generation mechanisms. 

The restriction to k > c is important as well. As proposition 0] shows, the 
behaviours of Efc(l/Z n ) and 1/M n can be quite different if k is not large 
enough. Proposition 0] deals with one generation of a branching process 
using random variables distributed as L c>mi when m — > oo, and corollary [21 
applies this result to the nth generation of a branching process using random 
variables distributed as L C)Tn for a given m, when n — > oo. 

Definition 3 Let Z denote a random variable and P^' m a probability mea- 
sure, such that Z is distributed, with respect to F^ m , like the sum of k i.i.d. 
copies of the random variable L c>m . 

Proposition 4 For any k ^ c, mK^ m (l/Z) — ► oo when m — ► oo. 
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Corollary 5 Assume that c is an integer and that £ n is the distribution of 
L c ,m, for every n. Hence M n = Zq m n . Then the distribution of Z n coincides 
with the distribution of the first generation of the branching process based on 
L c ,m n - As a consequence, m n E^' m (l/Z n ) — ► oo when n — > oo. 

Thus Efc(l/Z„) when k ^ c may decay on a different scale than 1/M n , see 
more on this case in section |HJ On the other hand, theorem [B] describes 
every square integrable Galton- Watson process if k is large enough, as the 
following theorem shows. 

Theorem C Any square integrable distribution on [1, +00 [ belongs to A c for 
c large enough, respectively to A' c for c small enough. Conversely, any distri- 
bution on [l,+oo[ which belongs to A c is square integrable and its variance 
is at most cm{m — \), where m denotes its mean. Likewise, the variance of 
any distribution on [1, +oo[ which belongs to A' c is, either finite and at least 
cm (m — 1), or infinite. 

We shall precise the optimal values of c for some usual distributions. Fi- 
nally, theorem |Bl indeed describes the behaviour of Efc(l/Z n ) when k is large 
enough, for any square integrable branching process. 

1.2 Bernoulli case 

We apply theorem [B] to the Bernoulli case when the offspring is always 1 
or 2. This case is relevant in the context of the polymerase chain reaction. 
Our techniques yield accurate bounds of Efc(l/Z n ) for every positive k, that 
is, even when k = 1, for instance in the homogeneous case, see theorem lEl 
below. We first state uniform bounds that are simple consequences of the 
results of section 11.11 

Theorem D Let n ^ 0. Assume that £j = (1 — Xj) 5\ + Xi 62 with Xi in [0, 1] 
for every i ^ n. Then, for any k ^ 2, 

1/k ^ M n E k (l/Z n ) ^ l/(k-l). 

In the homogeneous case, one can prove better bounds. We write E^ for E/% 
when £j = (1 — x) 5i + x52 for every i ^ 1. For every x in (0, 1), define 

a"(x):=l-x, a'(x) := (1 - x)/{l +x). 

Then ^ a' ^ a" ^ 1 and a" and a' decrease from a"(0) = a'(0) = 1 to 
a"(l) = a'(l) = 0. 
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Theorem E (1) For any k ^ 1 and n ^ 



l/k^(l + x) n EUl/Z n ) ^ l/(k 



a"(x)). 



(2) For any k ^ 1 



lim (l + x) n E£(l/Z n ) > \/{k 



a!{x)). 
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These estimates are precise enough to imply the following side result about 
the case k = 1. 

Proposition 6 There exists no uniform upper bound of (1 + x) n Ef(l/Z n ) 
for n ^ and x G (0, 1), since lim^o lim n _ >00 (l + x) n Ef(l/Z n ) is infinite. 
More precisely, for every x in (0, 1) and n ^ 1, 



where c(x) := 1 — x(l — x)/(l + 3x) is swc/i i/iat 8/9 ^ c(x) < 1. 

In theorem El the value of a'(x) stems from the general construction of sec- 
tion ll.ll but the value of a"(x) does not. In other words, a direct application 
of section n~T1 to the Bernoulli case yields a(x) instead of a"(x), with 



c(x)/x< lim (l + x) n Ef(l/Z n ) < l/x 



a(x) : 



log(l + x)/log(l — x). 



For every x in (0, 1), a'{x) < a"(x) < a(x). 
Theorem IE1 follows from the more general case below. 



Theorem F Let £j = (1 — Xi) 5\ + X{ 82 for every i. 



(1) If Xi ^ x for every i ^n, then, for any k ^ 1, 



n 



1/k ^ K k (l/Z n ) + Xi ) ^ l/(k - a"{x)). 




n 



lim E k (l/Z n ) Wl + Xi) >l/(k-a'(x)). 
1^00 x x 
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1.3 A discontinuity result 



In the Bernoulli case, the functions a"{x) and a'(x) have a nonzero limit at 
x — > + , hence the second part of theorem El above shows that the limit of 
the normalized harmonic moments does not always depend continuously on 
the parameters of the model. We show in this section that the phenomenon 
is general. For the sake of simplicity, we deal with the homogeneous case. 

Let M. denote a given subset of (1, +00) such that 1 is a limit point of M.. 
Below, the limits when fj, — ► 1 are implicitly restricted to m G Ai. For each 
M G Ai, let £ M denote a distribution of mean fj,. If & = £ M for every i ^ 1, 
define a function % on .M by 

Mm) := lim M"E fc (i/z n ). 



Proposition 7 Assume that, for each m € -M, i/iere exists o(m) and a '(M) 
siic/i i/iai £ M belongs to A a ^ and to TTien 

— a*) ^ liminf %(m) ^ lim sup ^(m) ^ l/(fc — a*), 

where a* := limsup^^ a(^) and a'* := limmf^_+i o'(m)- 77ms, i/ a* is 
positive, the function is not continuous at fj, = 1 + . 

Theorem G In £/ie homogeneous case, assume that each ^ is the law of 
1 + X , where the law of X is either binomial or Poisson or geometric. Then, 
for every k 1, is discontinuous at \x = 1, since fafc(l) = 1/k and 

lim h k (ji) = l/(k - 1). 
/x->l 

7f i/ie Zaiu 0/ X is geometric, then hk{n) = — 1) /or every \i > 1 and 

Likewise, assume that each ^ is the law of L c ^ for a given positive integer 
c. Then, /ifc(l) = 1/k for every k. Furthermore, for every fj, > 1, 7i&(m) = 
l/(k — c) if k > c and /ifc(/u) = +00 if k ^ a 



1.4 Higher harmonic moments 

We now state an extension of theorem iBl to higher harmonic moments. The- 
orem El is but a special case of proposition |^ in section [5J 

Theorem H (1) Let n ^ 1. Assume that there exists c such that £j G *4. c 
/or every i ^ n. Then, for every positive integer r and every integer k > rc, 

M; E fc (l/Z;) < l/[(fc - c)(fc - 2c) • • • (fc - rc)]. 
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(2) Assume that M n — > oo w/ien n — > oo and that there exists c such that 
£i G »4c /or every z. T/ien, for every positive integer r and every integer 
k > rc, 

lim A£E fc (l/Z£) > l/[(jfe-c)(ife-2c)---(A;-rc)]. 

n— »oo 

Corollary 8 Zei n ^ 1. Assume that there exists c such that & G »4 C /or 
every i ^ n and write af.{l/Z n ) for the variance of 1/Z n when Zq = k. 
Then, for every integer k > 2c, 

M 2 n a 2 k (l/Z n ) < (3c)/[fc(fc-c)(fc-2c)]. 

If, additionally, there exists d such that £j G ^/ /or every i ^ n, then the 
sequence 

k 3 M^a 2 k (l/Z n ) 

is bounded above and below by finite positive constants, independently of n 
and k, for large enough values of k. 

1.5 Related studies 

As mentioned above, Piau (2004) uses preliminary versions of our results, 
especially in the Bernoulli case, which is relevant for the study of the poly- 
merase chain reaction. In this specific case, we are now able to deal directly 
with every initial population k, even k = 1. 

Ney and Vidyashankar (2003) give asymptotics of the harmonic moments 
of every integrable homogeneous branching process starting from k = 1 
particle. When furthermore LlogL is integrable, their results specialize as 
follows, see also Bingham (1988) for some classical facts that are recalled 
below. 

Let p\ := P(L = 1), fi := E(L), and let 7 denote the Karlin-McGregor 
exponent of the distribution of L (7 is also called the Schroder constant), 
defined by the equality 

Pi /" 7 = 1- 

Let W denote the almost sure limit of the nonnegative martingale Z n //i n . 
The Poincare function is the Laplace transform P(s) := Ei(exp(— sW)) of 
the distribution of W when k = 1, and solves Poincare's functional equation 

p(H = f(p(s)). 

Three cases may arise. First, when r > 7, Ei(l/Z^)/p™ converges to a finite 
positive limit, whose expression is an integral which involves the Schroder 
function S, defined for any t in [0, 1), by 

S(t) := lim Ei(t z ™)/^. 
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Up to a multiplicative constant, S is the unique finite solution on [0, 1) of 
Schroder's functional equation 

S(f(t))= Pl S(t). 

Second, when r = 7, Ei(l/Zn)/(np") converges to a finite positive limit, 
whose expression involves Poincare and Schroder functions. Third, when r < 
7, fi nr Ei(l / Z^) converges to a finite positive limit. Ney and Vidyashankar 
provide an expression of the limit in terms of an integral of Poincare function. 
One can readily check that this limit is in fact Ei(l/W) and that the limit 
is also an upper bound. 

When L log L is not integrable, the results are similar but one must replace 
the normalizations rip™ = n/// 7 when r = 7 and 1/ ji nr when r < 7, by 
similar expressions which involve the Seneta-Heyde constants. 

Coming back definitely to the LlogL case, we recall that the distribution 
of W has a density w on the nonnegative real line, such that w(x)/x' Y ~ 1 
is bounded between two finite positive constants, when x — ► 0, see Dubuc 
(1971). 

The comparison of our results with those recalled above is based on two 
elementary lemmas. 

Lemma 9 For any distribution £ in A c , 7(£) ^ 1/c. For any distribution 
e in A' c , 7 (£) < 1/c 

In other words (see definition 221 in section 

a'(0 < 1/7(6 < MO- 

Lemma 10 For any branching process, k ^ 1, n ^ 0, and positive r , 

Corollary ^2 is not stated as such in the papers that we mentioned above 
but it follows from results that we recalled. 

Corollary 11 For any homogeneous branching process of Schroder expo- 
nent 7 and any k > r/7, the sequence ^ nr Efc(l / 'Z£) is bounded as n varies, 
by the finite constant Mi t (l/W r ). 

An interesting feature of corollary ^2 is that it deals with the entire regime 
where such a control of n n7 "Ek(l/ Z£) may hold, namely, with every initial 
population k > r/j(^). In other words, when k ^ r H{£)-> /^^(V-^n) * s 
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not bounded. Our upper bounds are restricted to higher values of k, namely, 
to the regime k > rA(£). 

One could think of recovering the dependence with respect to k from the 
results of Ney and Vidyashankar even when k ^ 2, starting from the in- 
equality 

E k (l/W r ) ^E 1 (l/W r/k ) k /k r . 

However, the bounds one gets cannot be optimal for k ^ 2, since 

E c k (l/W) = l/(k -c)< E\{l/W l,k ) k /k. 

Furthermore, as stressed by Bingham (1988), the law of W, hence the value 
of E(l/W r ), may be explicitly computed only in very specific cases. In 
contrast with every other paper we are aware of, the bounds we provide 
are explicit. The assumptions involve only elementary, step- by-step, mech- 
anisms of the branching process, that is, the distributions of the number 
of descendants at each generation. Also, we allow for inhomogeneous pro- 
cesses, as long as the reproducing laws belong uniformly to a given space 
A c , respectively A' c , and we make explicit the dependence of the bounds on 
the initial population. 

The introduction of the family of distributions described by g Ct m for integer 
values of c is hardly new, see Harris (1948) for instance. A key point is 
that we use them for noninteger values of c and as a reference scale of any 
square integrable distribution. For instance, the k = 1 case of Bernoulli 
distributions requires to make use of values of c in (0,1). Although these 
distributions do not correspond to a branching process for noninteger values 
of c, they still satisfy a semigroup property, and this property is sufficient 
for our purposes. Finally, our methods do not determine the behaviour of 
the harmonic moments of homogeneous processes whose reproducing law is 
not square integrable. 

1.6 Plan 

The remainder of the paper is organised as follows. In section [2J we reduce 
the case of branching processes in A c and A' c to the case of well-chosen 
distributions g c ,m, which we solve in section |31 In section 0J we show that 
the A c and A' c cases imply the result for every square integrable branching 
process. In section |S1 we deal with harmonic moments of higher degrees. 
In section El we study thoroughly the Bernoulli case, that is, the case when 
the offspring is 1 or 2, sharpening our previous results on this subject. We 
provide an algorithm to compute the asymptotic harmonic moments, up to 
any accuracy, and we present some simulations and conjectures about this 
specific case. Section [7] is a remark about size-biased offsprings. Finally, in 
section EJ we briefly explain how to deal with cases when the asymptotic 
behaviours of the harmonic mean and the mean do not coincide. 
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2 From g c m to A c and A' c 

We show that every branching process whose branching mechanism uses 
only laws in A c can be reduced to the case of L CyTn for a suitable m, and we 
solve this case. Similar results hold, as regards the comparison with A' c . 

2.1 Results 

Lemma 1121 describes the semi-group structure of each family (g c ,m)m- 

This 

is the starting point of our computations. Corollary is a special case of 
corollary and corollary E3 is a consequence of lemma 1121 Corollary ED 
uses definition 01 in the introduction. 

Lemma 12 For any positive c and any m ^ 1 and m' ^ 1, g c ,m ° <?c,m' = 
9c,m" with m" := mm' . 

Corollary 13 Let (p denote a nonnegative completely monotone function. 
For every branching process in A c , every k ^ 1 and n ^ 0, 

E k (<p{Z n )) ^ E% m {<p{Z)), where m := M n . 

For every branching process in A' c , every k ^ 1 and n ^ 0, 

E k ((p(Z n )) ^ E% m ((p(Z)), w/iere m := M n . 

Recall that <f> is completely monotone if and only if its derivatives are such 
that (—1)* is nonnegative for every positive integer i. Nonnegative com- 
pletely monotone functions are Laplace transforms of nonnegative measures 
on [0, +oo), see chapter IV of Widder (1948). 

Corollary 14 For every branching process in A c and every positive r, 
E k (l/Z r n ) ^E^' m (l/Z r ), where m := M n . 

For every branching process in A' c and every positive r, 

E k (l/Z r n ) ^ E£ m (l/Z r ), where m := M n . 

2.2 Proofs 

Proof of lemma 1121 Since each g C)Tn is the conjugate of g m by the bijection 
t i— > t c , the case c = 1 implies the general case. When c = 1, 1/(1 — g m (t)) 
is an affine function of 1/(1 — t). By composition, g m o g m , is also an affine 
function of 1/(1 — t) and it remains to compute its coefficients to prove the 
semigroup property. □ 
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Proof of corollary 1131 The representation of completely monotone func- 
tions which we recalled after the statement of the corollary shows that 

<p(z) = [ t z dn(t), 
Jo 

for a given measure n on [0, 1]. Thus E^(y(2 n )) is a positive linear functional 
of the generating function E k (t Zn ) of Z n . The function E k (t Zn ) is the fcth 
power of the composition from i = 1 to n of the generating function of £j. 
When the branching process belongs to A c , the generating function of £j 
is bounded above by gc,^-, thus the composition is bounded above by the 
composition of the functions g c ,/^> which equals g c , m - Finally, 

E k (<p(Z n )) ^ f 1 g Ctm (t) k dn(t) =E c l : m ^(Z)). 
Jo 

The proof of the result for branching processes in A' c is similar. □ 

Proof of corollary 1141 For every positive r, (p(z) := l/z r is completely 
monotone. To see this, choose div(t) = (log l/t) r ~ 1 dt/(T(r) t) in the repre- 
sentation of ip which we used to prove corollary 1131 □ 

3 The case g C)Tn 

Our task in this section is to evaluate the moments of 1 jZ under the measure 
P£' m . The cases k > c and k ^ c yield different asymptotic behaviours of the 
first moment of 1/Z. We begin with the direct way to deal with p^ ,m when k 
is large enough, namely, the computation of factorial moments of Z instead 
of the usual moments, see proposition 1151 Starting with lemma 1181 which 
gives a representation formula valid for every k, we study in depth the first 
harmonic moment, both in the small k and large k regimes. Corollary 1141 
in section |^1 and lemma below then imply theorem [Bj Lemma I^U] deals 
with the case k = c. Lemma I2T1 provides an alternative formulation of the 
integral of lemma 1181 a formulation that lemma uses to settle the case 
k < c. Proposition 01 which concludes the case k c, is then an easy 
consequence. 

3.1 Results 

We begin with exact formulas. Theorem [H] in section 11.41 is a consequence 
of corollary 1161 below. 
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Proposition 15 (1) For every nonnegative integer r, 

¥% m {Z{Z + c) • • • (Z + rc)) = m r+1 k(k + c) • • • (k + rc). 

(2) For every real number r such that k > rc, 

rrf¥^ m (T{(Z/c) - r)/T{Z/c)) = T((k/c) - r)/T{k/c). 

(3) For instance, for every nonnegative integer r such that k > rc, 

m r E jV*(i/[( Z _ c ) . . . ( Z _ rc )]) = !/p _ c ) . . . ( fc _ rc )j. 

Corollary 16 For every nonnegative integer r such that k > rc, 

l/k r < m r E^' m (l/Z r ) < l/[(k - c) ■ ■ ■ (k - rc)}. 
For instance, for every k > c, 

1/k ^ mE c /"(l/Z) < mE^' m (l/(Z - c)) = l/{k - c). 

Here is a slight generalization of the r = 1 assertion in corollary 1161 

Proposition 17 For every c>0,m^l,u^0 and positive integer k > u, 
1/k ^ mE^' m (l/(Z - «)) < l/(jfe - sup{c,n}). 

Proposition El and corollary El are the results that we use to settle the 
case k > c in the rest of the paper. We turn to the evaluation of the exact 
harmonic moment of Z with respect to P£' m . The results below are mostly 
used to deal with the case k ^ c. 

Lemma 18 For every positive integer k, positive c and m > 1, 

E£ m (l/Z) = G(k/c,m)/c, G(u,m) := f t^ 1 dt/(l + (m - l)t). 

Jo 

Alternatively, 

G{u,m) = B„,i_ u (l - l/m)/(m - 1)", 

where B U)V denotes the incomplete Beta function of parameters u andv, that 
is, for every x G [0, 1), 

b u>v (x) ■.= f e-^i-ty-Ut. 

Jo 
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Lemma 19 Assume that u > 1. Then, 

m G(u, m) ^ l/(u— 1). 

The order of this upper bound is exact when m is large, since the function 
(m — 1) G(u,m) increases when m increases and converges to l/(u— 1) when 
m — ► oo. 

Lemma 20 G(l, m) = (log m)/(m — 1). 

Lemma 21 For any u, (m— 1)" G(u,m) is an increasing function ofm ^ 1. 

Lemma 22 Assume thatu < 1. Whenm — ► oo, (m—l) u G(u,m) converges 
to c u := ir I sin(7ru). Thus, for any m > 1, 

(m - l) u G(u,m) < c u . 

On the other hand, for any m ^ 2, (m — l) u G(u,m) ^ l/(2u). 

Bounds of c u are, for every u < 1, c u ^ ir and l/(2n(l — u)) ^ c u ^ 

\/(u{l-u)). 

Corollary 23 Let k < c and £(k,c) := c/(k(c — k)). For every m, 

(m-l) fc / c E^' m (l/Z) <£(jfe,c). 
The order of this upper bound is exact, since 

lim (m- l) fc/c E'' m (l/Z) > £(jfe,c)/2. 

m^oo 

3.2 Proofs 

Lemmas I2UI 1211 and 1221 stem from the definitions. 
Proof of proposition IT5l (1) For any \x\ < 1 and any positive y, 
1/(1 -x)y = J2 x r T(y + r)/[T(y) T(r + 1)]. 

Setting y = Z/c and integrating yields 

E c fc < m (1/(1 - xf/ c ) = £ E* m [T(r + (Z/c))/T(Z/c)] x r /T(r + 1). 

On the other hand, 

E^ m (l/(1 - xf/ c ) = 5c , m (l/(l - x) 1 /^ = l/(l _ 

14 



Using the expansion of 1/(1 — mx) k l c given above and equating the coeffi- 
cients of the two series yield the result for any nonnegative integer r. 

(2) For any positive y and r with y > r, 

r(r) T(y - r)/T(y) = I t^' x (1 - if" 1 dt. 
Jo 

Setting y = Z/c and performing the integration yields 

r(r) n m \n(Z/c) - r)/T(Z/c)} = f g c , m (t 1/c ) k (1 - t)^ 1 dt/f+\ 

Jo 

The change of variables s := g c ,m{t l ^ c ) c = 9m (t) yields 

r(r) ¥% m [T{(Z/c) - r)/T(Z/c)} = f s fc / C -r-i (1 _ s) r-l ds/m r^ 

Jo 

that is, the desired formula. □ 



Proof of lemma HU Write E^' m (l/Z) as the integral of g c ,m(t) k /t on (0, 1). 
Use the change of variable t' := g c ,m(t) c - This yields the first expression of 
G in the lemma. To get the expression of G in terms of incomplete Beta 
function, use the change of variable t' := (m — 1)^/(1 + (m — l)t) in the first 
expression of G. □ 

Proof of lemma 1191 In the first expression of G in lemma 1181 use the fact 
that 1 + (m — l)t lies between rat and m. Thus, G(u,m) lies between the 
integral of t u ~ 2 /m and the integral of t u ~ l /m, that is, between l/((it — l)m) 
and l/(um). □ 

Proof of corollary 1231 Use the bound of lemma l22l bv c u /c for u := k/c, 

then the bound of c u by l/(u(l — u)). This yields the bound for every finite 
value of m. The limit when m — > oo is c u /c l/(2uc(l — u)) = £(k,c)/2. 
□ 



4 From A c and A' c to the general case 

In this section, we show that every square integrable branching process be- 
longs to the set A c , respectively to the set A' c , for a suitable value of c, we 
prove theorem Rland we describe the best possible constants c of theorem IE1 
in some specific examples. 
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4.1 Comparisons 



Our next proposition is related to theorem [0 and justifies definition 1251 
below. 

Proposition 24 Ifc\ C2 andm ^ 1, <7 C i,m ^ 9c 2 ,m- If ci < C2 andm > 1, 
the distribution of L C2)fn belongs to A C2 but not to A Cl and the distribution 
of L Cltm belongs to A' Cl but not to A 1 . Thus, (A c ) c is a strictly increasing 
sequence and (A' c ) c is a strictly decreasing sequence. 

Definition 25 For any square integrable distribution £ on [l,+oo[, let 
A(C) := inf{c >0;(G A c }, A'(£) : = sup{c > ; £ G A' c }. 

4.2 Examples 

We now study some specific transformations and examples. Proposition 1261 
follows from the definitions. 

Proposition 26 For every £, A'(£) ^ A(£). The inequality is strict except 
in two cases: either A(£) = -A'(£) = 0, and in that case, £ is a Dirac 
measure at m ^ 1; or A(£) = A'(£) = c is positive, and in that case, £ is 
the distribution of a random variable L Ct7n . 

Proposition 27 (1) If the laws of the independent 1 + X and 1 + X' belong 
to A' c , the law of 1 + X + X' belongs to A' c as well. This statement is false 
when one replaces A' c by A c . 

(2) If the law of 1 + X belongs to A c and if b is positive, the law of 1 + bX 
belongs to A c b- A similar statement holds if one replaces A c and A c b by A' c 
and A' cb . 

(3) If the law of L belongs to A c and if V dominates stochastically L, then 
the law of L' belongs to A c as well. For instance, if b is nonnegative, the 
law of L + b belongs to A c - A similar statement holds if one replaces A c by 
A' c . 

If £ is a Dirac measure, = = 0. Other usual cases are as follows. 

Proposition 28 (1) If £ is uniform on {l,...,n}, A(£) < 1. More pre- 
cisely, 2n A (€) =71 + 1. 

(2) //£ is uniform on {l,n}, A'(£) = (n - l)/(n + 1) and 2 A <^ +1 =n + l. 

(3) If £ is the law of 1 + X where X is binomial (n,x), A(^) < 1. More 
precisely, 

(l+xn)(l-x) nA ^ = 1. 
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(4) If £ is the law of 1 + X where X is Poisson of mean x, A(£) < 1. More 
precisely, 

With the notations of section IT31 cases (1) to (4) of proposition are such 
that Aft) = l/ 7 (6 > A'(0- 

To check that the three values A(£), 1/7(0 an d A'(£) can indeed be different, 
assume that f := (1 — p)5\ + (62 + 63) p/2 with p £ (0, 1). 

Then 7(6 := — log(l — p)/log(l + 3p). For p = 1/5, one can check that the 
function t 1— > K(t L )/gi^^ m (t) has positive derivatives at t = and t = 1. 
Thus some values of this function are greater than 1 and some are smaller 
than 1. This implies that A(£) > 1/7(6 > A '(0- 

4.3 Bernoulli case 

Definition 29 If £ = (1 — a;) 61 + x<52, write a(x) for A(£) and a'(x) for 
A'(0. 

Proposition 30 For any x £ (0, 1), a'{x) < a(x) < 1, since 

a'(x) = (1 - x)/(l + x), (1 - x) ai - x \l + x) = 1. 

Thus, a and a' decrease on (0, 1] /rom a(0 + ) = o/(0 + ) = 1 to a(l) = 
o/(l) = 0. -Bo£/i are discontinuous at since a(0) = c/(0) = 0. 

One can note that a'(x) ^ 1 — x ^ a{x). 

4.4 Proofs 

Proof of proposition 1241 Compare the logarithmic derivatives. □ 

Proof of theorem O Both results stem from the expansion of g c ^ m near 
1, which reads as follows, when t = o(l), 

5c,m(l - t) = 1 - mt + c (c + 1) m (m - 1) i 2 / 2 + °(* 2 )- 

On the other hand, 

E((l - = 1 - E(L) t + E(L(L - 1)) t 2 / 2 + o(t 2 ). 

A comparison of the second order terms of these expansions yield the con- 
dition on the variance of L for L to belong to A c , respectively to A' c . 
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To show that any square integrable distribution belongs to A c for suitable 
values of c, we first choose values of d and s < 1, both large enough to make 
sure that ~E(t L ) ^ gd,m(t) for every t ^ s. Thanks to the expansion above, 
this is possible for any d such that 

d(d+ l)m(m- 1) > E(L(L- 1)), m:=E(L). 

Then we choose a value of c > d large enough such that l/m 1 ' ^ gd.m(s)/s. 
Thus E(t L ) ^ gd,m(t) ^ g c ,m(t) for every t ^ s, and, since E(t L )/t is a 
nondecr easing function of i, for any t ^ s, 

E(t L ) ^ tE(/)/s ^ tg dtm (s)/s ^ t/m 1 ^ ^ g c>m (t). 

The proof for the comparison with distributions in A' c is similar. □ 

Proof of proposition 1271 Part (1) stems from the property 

gc,m(t) gc,m' it) t gc,mm' (0 ; 

which we leave as an exercice for the reader. Parts (2) and (3) are direct. 
□ 

5 Higher moments 

Assume that £j G A c for every i ^ n and let m := M n . Then, 

Mr n M k (l/zr n )^m r E c > m (l/Z r ). 

Expansions of g c ,m(t) when m —* oo show that the distribution of Z/m with 
respect to P^' m converges to the distribution of W with respect to a measure 
Pf such that 

El(e- tw ) = (l + ct)-^ c . 

The distribution of W is Gamma (c, 1/c), that is, its density with respect 
to the Lebesgue measure dw is 

Furthermore, g c , m {t) ^ g c {t) := Ef(t u/ ) = (1 - c logt) -1 / c . Hence, 

< ^.(V^ r ) = m/c) - r)/(c r F(k/c)). 

This inequality holds for every positive values of r and c and every positive 
integer k such that k > cr. The lines above prove the following result. 
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Proposition 31 (1) Let c such that £j 6 ,4 C for every i ^ n. For every 
positive r and k such that k > r c, 

M T n ¥, k {l/Z r n ) ^ l/[Jfe,c] P , w/iere [fe,c] P := c r T(k/c)/T({k/c) -r). 

When r is an integer, 

r 

[k,c] r = Y[(k - ic). 
1=1 

(2) Conversely, let c such that & £ ^ /or every i. Assume that M n — > oo. 
Then, for every k and every positive real number r ^ k/c, 

lim M;E fe (l/Z;) = +oo. 

n— >cx> 

6 Bernoulli branching processes 

6.1 Preliminaries 

We first set some notations, to be able to deal with non homogeneous pro- 
cesses. 

Definition 32 The efficiency of a Bernoulli branching process is the se- 
quence X := (xj)j^i such that = {1 — Xi) 5\+ Xi 82- Let L, respectively L* , 
denote the set of efficiencies X such that X{ G [0, 1], respectively Xi £ (0, 1], 
for every i ^ 1. For any X £ L, let s(X) := denote the shifted 

sequence. 

Definition 33 For any k ^ 1 and any efficiency X , let 

n 

B k {X) := lim E fc (l/Z n )TT(l + ^). 

n— >oo -*- -*- 

i=l 

Ln the homogeneous case Xi = x for every i ^ 1, we write B k {x) for B k (X). 

By convexity, the limit which defines B k (X) is also a supremum over n ^ 0, 
thus B k {X) ^ 1/k. The functional B^ describes Efc(l/Z n ) for finite values of 
n as well, since replacing every Xi with i ^ n + 1 by freezes the branching 
process at its value Z n . Thus, uniform upper bounds of B k on L yield upper 
bounds of E^.(l/Z n ) for finite values of n. 

6.2 Results 

The following uniform result is a consequence of the fact that A(£) < 1 for 
every Bernoulli £, see proposition 
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Proposition 34 For every efficiency X £ L and every k #S 1 



1/jfc < 5 fc (x) < l/(jfc - 1). 



27ms, /or every n ^ 0, 



1/k < E fc (l/Z n ) JJ(1 + x,) < l/(fc - 1). 



i=l 



The sequence (-Bfc)fc>i satisfies recursion relations, which we state in propo- 
sition and which characterize it fully, see proposition I3fil 

Proposition 35 For every k ^ 1, f/je function B k is measurable on L. 
Furthermore, for every k 1 and <Y £ L, 



Proposition 36 Let (-Pfc)fe>i denote a sequence of functionals defined on 
L* . Assume that k F k (X) — > 1 when k — ► oo, uniformly over X £ L* , and 
that (-Pfe)fe>i solves (>EJ) on £* /or every k 1. These conditions define a 
unique sequence (i*fc)fc^i, suc/i i/iai i*& = on £* /or every k 1. 

Thus, the sequence (-Bfc)fc>i is entirely determined on L* by the recursion 
and by the bounds 1/fc < B fc < l/(fc - 1). 

Note finally that the recursion @ is but a special case of the following result. 
For any branching process of reproducing law S = and any k ^ 1, 

introduce 



and the shifted mechanism s(H) := Let denote the convolution 

of the measure £i with itself k times. Then, 



6.3 Homogeneous case 

We start with a version of the relation (J*j) in the homogeneous case. 
Proposition 37 For every x in (0,1), 




(*) 



H k (E) := lim M n E k (l/Z n ) 
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The recursion whose left hand side is Bf.{x) involves the whole set of values 
Bk(x), Bf, + i(x), . . . , B 2 k(x). Thus, this system of equations does not yield 
directly the value of each Bk(x). The exception is the case k = 1. 

Corollary 38 For every x ^ 0, B\{x) = B 2 {x) (1 + x)/x. 

Our main result in this section is proposition 1391 

Proposition 39 Leta'(x) := (1 — x)/(l + x) anda"{x) := 1 — x. For every 
fe > 1, 

l/(fc - a\x)) ^ B k (x) < l/(k - a"(x)). 
Thus, B k (0+) = l/(k - 1) and B k (l~) = 1/Jfe. 

For ft = 1, propositionlBUlstates that B\{x) is at least 1/(1 — a' (x)). A better 
bound obtains if one uses corollary 1381 and then proposition \\\\)\ namely 

(l + x) 2 /(l + 3x) < xB t {x) < 1. 

The lower bound is always greater than 8/9 = .889". From our numerical 
simulations in section [6.91 below, some values of \B\(X) are as small as 
B* = .9274 ± .0002, to be compared to 8/9 = .8889". 

One could iterate the procedure, getting yet tighter upper and lower bounds 
of Mf(l/Z n ), or of any ¥%(1/Z n ) with k 1, with any prescribed accuracy. 
We develop this idea in section IfTTH below. 

We end this section with a conjecture. 

Problem 40 We conjecture that every function x ^ Bk{x) is decreasing on 
x £ (0, 1]. Prove this and find a natural explanation of the fact that B k (0 + ) 
and -Bfc(O) are not equal. 

6.4 Proofs 

Proof of proposition l35l Since each Kk(l/Z n ) is measurable with respect 
to (xi)i<^ n , By.[x) is the limit of a measurable nondecreasing sequence, hence 
Bk is measurable. As regards the recursion relation, we consider the con- 
ditioning by Z\ of the Bernoulli branching process starting from Zq = k. 
On the event {Z\ = k + i}, (Z n+ i) n ^o follows the law of a Galton- Watson 
branching process of efficiency s(x), starting from k + i. Hence Q follows 
from the fact that the distribution of Z\ — k is binomial (k,x±). □ 
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Proof of proposition 1361 The existence follows from the construction of 
each B^. A proof of the uniqueness is as follows. Assume that the sequences 
of functional (B' k ) and (B k ) are solutions. In particular, B' k /B k — > 1 uni- 
formly on L*, when k — > oo. Fix e. For every k large enough and for every 
x G L*, 

B' k {x)^{l + e)B' k \x). 

Since (B' k ) and (B k ) both solve the recursion relations a recursion over 
the decreasing values of k shows that B' k (x) ^ (1 + s) B' k '(x) for every k ^ 1 
and every x £ L*. This recursion uses as a crucial tool the fact that no Xk is 
zero. Now, since e is arbitrary, B' k ^ on L* for every k ^ 1. Exchanging 
the roles of the two sequences, one sees that B' k = B' k on L*, for every k ^ 1. 
□ 

6.5 Outline of the proof of proposition 1391 

We start from relations between the functions B k in proposition 1371 which 
read, for every k ^ 1, 

B fc (x) = (l + x)Ejg(S Zl (x)). 

With respect to the probability P?, Zi is distributed like the sum of k i.i.d. 
random variables of distribution (1 — x) 5% + x oV Lemma HTl below follows 
from the fact that kBk(x) — > 1 when k — > 00. 

Lemma 41 Assume that liminf kip(k) ^ 1 and i/iai, /or every k ^ 1, 

(1^)EJW))<#). 

T/ien Bk{x) ^ / or every k ^ 1. Conversely, if limsup /c ^ 1 anc ^ 
i/j /or every k ^ 1, 

then Bk{x) ^ V'(^) / or every k ^ 1. 

Definition 42 For every k ^ 1, let Ck denote the unique solution in (0, 1) 
0/ i/ie equation 

(l + x)E^(l/(Z 1 -c fc )) = l/(fc-c fe ). 

Lemma |43l becomes obvious when one uses an equivalent definition of Cfc, 
given below in part (ii) of lemma Hoi 

Lemma 43 // c& ^ c /or every k ^ 1, then B k (x) ^ — c) /or every 
k ^ 1. Conversely, if Ck ^ c /or every k ^ 1, then Bk(x) ^ l/(k — c) /or 
every k ^ 1. 
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Lemma H31 asserts that Bk(x) ^ l/(k — c) for c = supjcfc ; k ^ 1} and, by 
lemma IHl this supremum is c\ = 1 — x, thus proposition 1391 follows. 

Lemma 44 For every k 1, ^ c\ = 1 — x. 

The following result shows that the technique above cannot yield a better 
value of a'(x) than a'(x) = (1 — x)/(l + x). 

Lemma 45 When k —* oo, c& — ► (1 — x)/(l + x). Furthermore, for k ^ 2, 

Cfc + xc fc _i ^ 1 - x. 

Finally, we use the characterizations below to evaluate c&. 

Lemma 46 For every k ^ 2, i/ie following inequalities are equivalent and 
equivalent to the fact that c ^ c& . 

(i) (l + ^E^l/CZi-cJXl/Cfc-c). 

(ii) fc (1 + z) E^_ 1 (1/(Z 1 + 2 - c)) > 1. 

(iii) fc (fc - 1) x (1 + x) E|_ a (l/(Zi + 4 - c)) < x k - 1 + c. 

The reversed inequalities (i'), (ii') and (iii') are equivalent and equivalent 
to the fact that c ^ . 

6.6 Technical steps of the proof of proposition 1391 

We prove lemmas El an d 14 5 1 assuming lemma HH1 for the moment. By 
Jensen's inequality, the expectation of the inverse is greater than the inverse 
of the expectation. Thus (ii') implies 

k{l + x) ^ (k - 1)(1 + x) + 2 - c. 

This reads c 1 — x. Since c\ = 1 — x, we are done with lemma 1441 
Furthermore, we can and we will restrict the reasoning below to c ^ 1 — x. 

To prove lemma I45| we first note that (ii) involves the expected value of a 
concave function of u := l/(-Zi — Ck-i), namely the function u i— > u/(l + bu) 
with 6 := c&_i + 2 — c. The expected value of a concave function is at most 
the function of the expected value. From the definition of c^-i, (ii) implies 

fc(l + x) > (k - 1 - cfc_i)(l + s) + c fc _i + 2 - c. 

This is equivalent to c ^ 1 — x — x Ck-i- Hence, for any k ^ 2, 

1 - x - x c fc _i < c fc ^ 1 - x. (f) 
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This is enough to show that ^ (1 — x) 2 for every k 1. Thus we can and 
we will further restrict the reasoning below to c ^ (1 — x) 2 . 

In the second step of the proof of lemma H51 we use (hi') like we used (ii). 
Namely, we note that (hi') involves the expected value of a concave function 
of l/(Zi — Ck-2) and we apply Jensen's inequality once again. From the 
definition of (hi') implies 

k (k - 1) x (1 + x) ^ [(1 + x) (k - 2 - c fc _ 2 ) + c fc -2 + 4 - c] (xk - 1 + c). 

After some simplifications, this reads + Aq ^ 0, with 

^1 := (1 - x) 2 - c + x 2 c k -2, A := (1 -c)(2(l - x) - c - xc k - 2 )- 

Since c ^ (1 — x) 2 and Cfc_2 ^ (1 — x) 2 , simple bounds show that Aq ^ 1. 
Hence (hi') implies that Ai ^ —1/k. Finally, for every k 3, 

(l-x) 2 ^C fc ^(l-x) 2 +X 2 Cfc_2 + lA. (J) 

One uses the a priori bounds of (m) and (jB as follows. On the one hand, the 
upper bound of c k in ((|J) implies 

lim sup Cfc ^ (1 — a;) 2 + x 2 lim sup c&. 

On the other hand, the lower bound of in Q implies 

1 — x — x lim sup Cfc ^ lim inf . 

Hence limsupcfc = lim inf c k = (1 — x)/(l + x). This proves lemma H*5l 
Lemma HE] is a consequence of the following trick. Part (i) involves 

(k - c)/(Zi - c) = 1 - {Z x - k)/(Z 1 -c)=:l-v. 

By exchangeability, E£(f) is A; times the expected value of (Li — 1)/ (Zi — c), 
where Li denotes the number of descendants of the first individual in the 
initial population. The event {L\ — 1 7^ 0} is {L\ = 2} and has probability 
x. Thus, for every k ^ 2, 

(k - c) E£(1/(Z! - c)) = 1 - kxEUil/iZ, + 2 - c)). 

With the convention that Z\ = 0, P§ almost surely, this relation holds for 
k = 1 as well. This translates (i) or (i') into (ii) or (ii'). The translation of 
(ii) or (ii') into (hi) or (hi') uses the same trick, starting from 1/(Z\ + 2 — c). 
This concludes the proof of proposition EM 
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6.7 Algorithm 

The following algorithm yields approximate values of on L*, with any 
prescribed accuracy. 

• Fix n 1. 

• For every k n + 1 and x, let 

Sg (ft (x) := Bl n (x):=l/(k-l). 

• Find the unique sequence {B\ n )k^n that solves the system of equations 

for k ^ n, when one replaces every Bk(s{x)) such that fc ^ n + 1, 
by Bl n (s(x)), that is, by the value l/(fc — 1). 

• Likewise, find the unique sequence (B® n (x))/ c <g n that solves the system 
of equations (J*j) for k ^ n, when one replaces every Bk(s(x)) such that 

^ n+ 1, by B\ n (x), that is, by the value 1/k. 

• Then, for every k ^ 1 and every x, 

Bl n (x) ^ B k {x) ^ Bljx) < (1 + 1/n) fl° n (x). 

6.8 Comments on the algorithm 

Neither (-B^ n )fe>i nor {B\ n )k^\ solve the full system of equations For 
any fixed values of k and x, (-Bfc n (^))n^i is a nondecreasing sequence that 
converges to -Bfc(x) when n — > oo, Likewise, (-Bfc n (^))n^i is a nonincreasing 
sequence that converges to -Bfc(x) when n — > oo, 

Increasing values of n yield more and more accurate approximations of each 
-Bfc(x) and the relative error is of order at most 1/n. 

In the Bernoulli case, one can use some initial values, better than 1/k, 
respectively l/{k — 1), namely, for every k ^ n + 1 and x, 

b° k>n (x) := l/(k - a'(x)), b l Kn {x) := l/(k - a"{x)). 

The relative error that was at most 1 + 1/n in the first version of the algo- 
rithm becomes at most 

1 + (a" - a')(x)/(n + 1 - a"{x)) < 1 + (3 - 2\[2)/{n + x). 

Numerically, this is at most 1 + .172/n, for every x. 
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6.9 Simulations in the homogeneous case 

The algorithm above with n := 1000 suggest the following refinements. 
Define 

B(x) := x = B 2 (x) (1 + x). 

Simulations show that B decreases on (0, x*) from B(0 + ) = 1 to B{x if ) := B* 
and increases on (x*, 1] from B* to B(l) = 1, with 

x* = .38 ± .01, = .9274 ± .0002. 

This would imply that, for every positive x, 

B*/x ^ Bx[x) ^ 1/x. 

Simulations show that B 2 , hence B\ as well, decreases on (0, 1]. 

7 Size-biased offspring 

When computing harmonic means, it may prove convenient to use size- 
biased distributions, as follows. Assume that L and Lj are i.i.d. positive 
integrable random variables and that L' is an independent size-biased copy 
of L, that is, for every t G [0, 1], 

E(t L ') :=E(Lt L )/E(L). 

Then, for any nonnegative integer k, 

E(L) E(l/(Li + • • • + L k + L')) = l/(k + 1). 

Can one use this in our branching setting? Assume first that 1 ^ L ^ c + 1 
almost surely, for a given integer c. Since L' ^ c + 1 ^ L^+i + ■ ■ ■ + -£/c+ c +i 
almost surely, this proves that 

E(L)E fc (l/Zi) < l/(k-c), 

for every ^ c + 1. More generally, the last inequality above holds as soon 
as ~E(t L ) > E(t L ) c+1 for every t G [0, 1] and for a given positive c. 

In our setting, this line of reasoning suffers from two drawbacks. First, unless 
we miss something, to be able to iterate this inequality during n generations, 
one must assume that k ^ nc. Second, the inequality E(i L ) ^ E(i L ) c+1 
implies that c ^ -4(0 > where £ denotes the law of L (the proof is easy and 
omitted). In other words, k ^ c implies that £ G -4 C . 
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8 Case k < c 



This section is a brief description of the behaviour of Efc(l/Z n ) when the 
hypotheses of theorem iBl are not met. Consider, for the sake of simplicity, 
a homogeneous branching process and let 

Pi := £(*) = HL = %)= P(^i =i\Z = l). 

Our first remark is that, for every n ^ and k 1, 

E fc (l/Z n ) > r£/fc, r fc := max{rf, 1///}. 

The l//x bound is due to the convexity. The p\ bound is due to the fact 
that the probability of the event {Z n = k} is p\. 

The parameters \x and p±, through rjt, indeed describe the asymptotics of 
Efc(l/2/ n ), as follows. For the sake of simplicity, we exclude the degenerate 
case p\ n = 1, where polynomial corrections appear. For a given k ^ 1, 
there exists a finite positive hk such that 

lim E fc (l/Z n )/r£ = h k . 

In the Bernoulli case, p\ < 1/fi hence r& = for every k 1 and the 
r fc = Pi regime is nonexistent. 

The limits hk satisfy the following relations. Assume for instance that one 
wishes to compute hi and that fipi > 1, hence r± = p\. Conditioning on 
the value of Z\, one gets a relation between Ei(l/Z n+ i) and the sequence 
(Ek(l/Z n ))k^i- Letting n go to infinity yields 

h 1 = l + ^2H hPk /p 1 , H k :=Y J ^ka/Zn)M- 

The term ^(l/Z^/p™ behaves like (r^/pi)" - , that is, like l/(/ipi) n if \ip\ > 
1, or like if \xp\ < 1. Since both quantities are summable, H2 is finite. 
Since Hk ^ H2 for every k ^ 2, h\ is finite as well. 

When [ip\ < 1, = for every k ^ 1 and the same reasoning as above 
yields 

hi = H ^Pk hk- 

Since /ipi < 1, this gives hi as a linear combination of {hk)k>2- In turn, for 
every ^ 2, a one-step recursion, similar to the one we used before, shows 
that hk is such that 
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Since the coefficient of hk in the right hand side is [ip\ < 1, this relation 
implies that hk is a linear combination of the sequence (/ii)i>fc+i and that 
this linear combination uses nonnegative coefficients. Since hi ^ /12 for 
every i ^ k + 1, this series converges. However, it does not seem easy to get 
information about the coefficients h\ or hk with k ^ 2 from these relations. 
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